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Abstract 

In this article, we have introduced a VT symmetric non-Hermitian Hamihonian model which is 
given as H = uj{Vh + i) + a(6^ - (6''")^) where cj and a are real constants, b and tit are first order 
differential operators. The Hermitian form of the Hamiltonian T-L is obtained by suitable mappings 
and it is interrelated to the time independent one dimensional Dirac equation in the presence of 
position dependent mass. Then, Dirac equation is reduced to a Schrodinger-like equation and two 
new complex non-VT symmetric vector potentials are generated. We have obtained real spectrum 
for these new complex vector potentials using shape invariance method. We have searched the real 
energy values using numerical methods for the specific values of the parameters. 
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1 Introduction 

The nature of quantization arises due to the symmetry of the equations governing the physics. So, 
symmetry has long been known as a powerful and computational topic in quantum mechanics. Recently, 
there have been many studies of VT symmetric non-Hermitian systems with real energy since the original 
work of Bender and Boettcher [T] and the literature on such systems has expanded rapidly [11131 HI [SJ [SI [Z] • 
One of the key points of the investigation was the VT symmetry generated by the product of the 
parity, V, and time,7~ , linear and anti-linear inversion operators VxV — —x, TxT — x, TiT — —i. 
The operator T is anti- linear because it changes the sign of i. If VT symmetry of the Hamiltonian 
is unbroken; eigenfunction of the operator VT is simultaneously an eigenstate of Hamiltonian H , i.e. 
[H,VT] = 0. Later, it has been realized that the existence of real eigenvalues can be associated with 
a non-Hermitian Hamiltonian provided it is ?7-pseudo- Hermitian '5': r/H — Hh] where ?7 is a Hermitian 
linear automorphism which can be given as ?7 = (OO^)^^, O is a linear invertible operator. Here, the 
Hilbert space equipped with the inner product < .,77. > is identified as the physical Hilbert space. And 
the observable & which is the element of physical Hilbert space is related to the Hermitian operator 6 
by means of a similarity transformation Q — p~^9p where p = rj^ . At the same time Bagchi and Quesne 
have established that the twin concepts of pseudo-Hermiticity and weak-pseudo-Hermiticity [5]. Thus, 
the concept of pseudo-Hermiticity has attracted much interest on behalf of physicists [lOl [11] [121 [S] • 



^e-mail : yesiltasOgazi. odu.tr 



1 



In T?', the Kepler problem solutions are investigated in Dirac theory for the particle whose mass 
is position dependent and the effective mass is given in the form of a multipole expansion, existence 
of the bound states are discussed in detail. Earlier, using a standard expansion of radial functions as 
a different approach, Dirac equation spectrum was obtained for the mixed potentials [TS]. Results of 
[14) which are about large quantum numbers not leading to inverse mass and momentum independent 
energy is also consistent with those found in [TB]. In [T^, Dirac matrices ci,/3 have space factors as 
/, /i functions where / responsible of the deformation of the Heisenberg algebra for the coordinates 
and momentum operators, /i is responsible of a dependence of the particle mass on its position. Exact 
solutions are found for the fermion in Coulomb field with the function / which depends on r linearly while 
the function /i depends on r inversely. It is pointed out that the spectrum results of [17] can be useful 
for the nanoheterosystems. Similar arguments about the Dirac oscillator with deformed commutation 
relations leading to the existence of the minimal length of space can be found in jTS) . 

Latterly, non-Hermitian potentials for the fermions have been studied in the literature [THl dHl HIl EH 
123] and fermion models interacting with VT symmetric potentials in presence of effective mass have been 
attracted interest [24l [25l [26l [27l [28l [29l [30] . In [31] , the one-dimensional effective mass Dirac equation 
bound states are studied within the interactions of non-PT-symmetric, and non-Hermitian, exponential 
type potentials. Moreover, (1 -I- 1) Dirac equation with position-dependent mass (PDM) and complexified 
Lorentz scalar interactions, is discussed through the supersymmetric quantum mechanics |32| . More 
references about the complex potentials and Dirac theory can be found in [3^ . 

Also, pseudo-Hermitian interaction in relativistic quantum mechanics is studied with the positive 
definite metric operator rj calculations for the state vectors [34]. Using the spin and pseudo-spin concept, 
spectrum of PT symmetric Rosen-Morse potential is studied and analytical methods are used in |35j . 
Dirac equation with position dependent effective mass transformed into Schrodinger-like equation is 
studied in a general context and Levai's method is used |36| . Supersymmetric quantum mechanics 
(SUSY QM) provides elegant procedures to solve some classes of potentials with unbroken SUSY and 
shape invariance (SI) which is one of the standard way and it is known that the potential algebras of 
these systems have been investigated to find exact solutions [371[3Hl[MllinillIlllllMllllll3llS]- SUSY 
QM methods and relativistic extensions have been used by many authors [IT] |3H1 HSl [5^ ■ 

The purpose of the present paper is to explore new relativistic complex vector potentials of the non- 
Hermitian bosonic Hamiltonians which may be unsolvable and map them into a solvable but real effective 
potentials. In the literature, bosonic/fermionic Hamiltonians with two mode have physical importance 
such as Jaynes-Cummings model in solid-state physics [51], Bose-einstein condensate [52], squeezed states 
in a condensate of ultracold bosonic atoms confined by a double- well potential [53] . 

Using the methods of SUSY QM, we have obtained solutions of complex vector potentials and showed 
that in Dirac equation, decomposing the the vector potential into the real and imaginary parts leads to 
derive both exactly and conditionally exactly solvable potentials. The paper is organized as follows: In 
section 2 a non-Hermitian Hamiltonian model is introduced by us and mapped into its Hermitian form. 
Shape invariance which is one of the effective tool in SUSY QM is given shortly in section 3. Section 4 
includes the mapping of the Dirac equation into a Schrodinger-like equation and obtaining new complex 
and effective potentials with their exact solutions. 

2 The Non-Hermitian Model and Hermitian Equivalents 

Previous works by the authors have included many aspects of a non-Hermitian su{2) Hamiltonian known 
as Swanson Hamiltonian [T^ [M] [SS] [SHI [ST] [55]. The Swanson Hamiltonian is given by H = aj(a^a + 
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1/2) + aa^ + j3a^ , where a, at are annihilation and creation operators, w, a and /3 are real constants. In 
this paper, let us consider a VT symmetric non-Hermitian model with two parameters given by 

n = Uj{Pb+]^)+a{p -(Pf) (1) 

where f is Hermitian adjoint, h is the annihilation operator given in a general form, 

h = A{x)^+B{x) (2) 
ax 

and A{x), B{x) are real functions. VT operator has an effect as a; — > —x, p ^ p and i — — i in the 
Hamiltonian, if the operators are taken as 6 = ^a; + -^P, it can be seen that the Hamiltonian is VT 
symmetric. Now, in terms of differential operators, ([Ij becomes 

H = ~ujA{xf— + UaA(x)B(x) - 2ujA(x)A(x) ) — 

ax^ ax (^3^ 

-{lu - 2a)A{x)B{x)' - (cj - 2a)A{x)' B{x) + ujB{xf - a(A(x)A(a;)" + {A{x)' f) + |. 
We may write the eigenvalue equation for ([TJ as given below 

-Hi! = eV- (4) 

Here, the pseudo-Hermitian Hamiltonian ([3]) can be mapped into a Hermitian operator form by using a 
mapping function p 

h = pUp-^ (5) 

where 

Here we note that hifj = £-0, jjj = p~^^. So we can introduce operator h which is Hermitian equivalent of 
v. as 

^ = -''Tx^^''^'Tx^^'f'^'''^ 

here Ueff{x) takes the form 

U^fjix) = I - u;iAix)B{x))' a [{A (x))' + A{x)A" (x)) + + ^)B'{x) (8) 
where the primes denote the derivatives. Then ([7]) can be mapped into a Schrodinger-like form by using 

A{x) 

Hence, Schrodinger-like equation becomes 

-*'V) + fe-^«2: + ^|M + — 4M_;^(4^ ,10) 

\iuA^{x) A^{x) uj^ A^{x) LU A{x) uj A^{x) 
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3 Shape Invar iance 



It is very well-known that a quantum system having a square- integrable ground-state with finite/infinite 
discrete energy levels Eq < Ei < E2 < ... where the ground-state energy is chosen to be zero £^0=0 
is a fundamental idea in supersymmetric quantum mechanics. Generally we can denote the positive 
semi-definite Hamiltonian by H can be given in a factorized form |46] : 

m = AU = —^ + v{x) (II) 

A=^~W{x), A^ ^-^ + W{x) (12) 
ax ax 

v^{x) ^W^{x)±w'{x). (13) 

We used the unit system h = 2to = I. Here W{x) is the function which is real and smooth known as 
the superpotential and the ground-state wave- function Co (a;) = e" '^^^^f^/) ig nodeles. It is noted that 
ACoix) = 0. In this approach, potential depends on a set of parameters a — (ao, ai, 02, ...) to be expressed 
by W{x, a),A{a), E{a), ... The shape invariance condition is 

A{a)A\a) = A^a + A)Aia + A) + Ei{a) (14) 

in which A is the shift of the parameters. The entire set of discrete eigenvalues and corresponding 
eigenfunctions are En{a) and Cn{x, a) can be written as 

n-l 

En{a) = ^£;i(a-|-fcA) (15) 

fc=0 

C„(x,a) - A\a)A'^{a + l^)...A\a+{n-l)l^)e~^^'^y^^y^''+''^\ (16) 



4 Dirac Equation 

The Dirac equation which plays an important role in relativistic quantum mechanics describes relativistic 
effects due to the speed and spin of particles. The one dimensional time independent Dirac equation with 
effective mass M[x) and vector potential V{x) is 

{a.p + j3M{x) + Vi)'^{x) = eN{x) (17) 

where ^ is the two component spinor wave-function, E is the energy, pis the momentum operator, M{x) 
denotes the position dependent mass and a and (3 are 2x2 Dirac matrices in standard representation and 
h = c = I atomic units are chosen. Let us show the upper and lower components by 4'{x) and 6{x). Using 
a = (73, /? = tJi, where cti and are Pauli matrices, and multiplying (jI7l) by cri, then we obtain |I9| 

+ (E -V(x))e - M(x)(j) = 

ax 

i^ + (E-V(x))^-M(x)e = 0. 
ax 
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If we terminate 9 in above coupled differential equations, we obtain 

(18) 

We use a transformation of the upper component wave- function which is 4>{x) = M{x)ip{x) in (|18p . we 
find that ^ 

+ = (19) 

Here, effective potential Veff{x) reads 



^e//(a:) = -V^{x) - z—^ + M\x) + ^-y-— -Li + E 2V{x) - ^ > 



dx M{x) dx V ^^{^) 

1 d2jvf(2;) 3 / 1 dM{xy " 



(20) 



2M{x) dx^ 4 VM(x) dx 
Now we decompose the vector potential V{x) into the real and imaginary parts in (|20p as 

l^(a:) = V^i^(x)+zF/(x) (21) 

which leads to 

9/ X 9/ X 9/ X /X M" (x) 3 f M'(x)\^ ^ M'(x) , , 

V.^x) = -V-(.) + Vhx) + Af^(.) + 2i.l...(^) - ^ + i j + (-) - -J^Vd^) 

(^-2V^/(x)l^;^,(^) + -^EViix) ~ F^(:.) + ^^Vi,(x) - E^^^ . 

(22) 

We may terminate the imaginary part of V^f / (x) by using 

^ M(xr ^h(^) .„x 

^ 2Af(x) 2(£; - ^flXa;)) ■ 

Because we have obtained a real effective potential expression for the non-Hermitian Hamiltonian in the 
last section. Now, we can give Veff{x) in the form of 

In order to compare Ueff{x) and Ve//(x), we may choose M(x) and VR(a;) as 

MW = + (25) 
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and put in (|24p where mi and m2 are real constants. Thus, we give another ansatze for B{x) as 

B{x) = jA{x)+PA'{x) 
where 7 and (3 are real constants. Afterwards, Veff{x) takes the form given below: 

, 22, o.,™ , «™ ^^(^)' , fr.^ , «^ ^2 l^ f Mx)' Y , 



This time, we shall use pT]) in ([T0| so that wc would compare (|28l) and (ITO| . then we obtain 



a;2+4a2 ^ uj/2-e ( ^u? + a\{A(x)f 



- /? + 27 /3 - 1 = e$(a;) 



w / A{x) \ Lo'^ J A{x) 

and we can also give Uef / [x) as 

Ueff{x)^ r+e+ \, .0 +\P 5 P- 



A{xY \ w2 ^ ) A{xY 



A{x) \ uj"^ ) A{x) ' 

Hence, we can compare and (l30l) and ((28|) . then we find this set of equations 

22 + 2 
e = 7 TO2 2 7 



2w 
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ni2[mi + pm2j — 5 P ~ 1- 



From the last relation we can find 
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nil = — {~Pwm2 ± i/a;2(l + (S^m^) — 4auj) 
and then, we can give E in terms of parameters oj, a as 

2 2/^2 + 

^ = 2~^ r^^^;^ 

Now we will give two potential models: 



4.1 Example 1: non-PT symmetric vector potential 

Using some special values of A{x) may give rise to solvable effective potential models. For instance, if 
A{x) — 5 cosh a; is chosen, one obtains 

V^(x) = £ - £ sec /ix + — — (37) 

2 /i cosh X + sm hx 

that is not a solvable non-VT symmetric potential, at the same time, the mass expression is given by 

Af(x) = TO27 + TO^^ l3~l^ta.nhx. (38) 

In this case Veff{x) is obtained as 



Veff{x) ^ - yE'^ ^ 4 ^ '•"^^ +/3w2) j sec/ix +27TO2(»7ii+/3"^2)tanha; + 7'^m2 + - + (mi+/3m2) . 

(39) 

We can give (1591) in terms of w and a constants by the aid of (PT|) - ([M1) : 

Vef f {x) = Vb — Vi sec h^x + V2 tan hx, — 00 < x < 00 (40) 

where 

uj , 1 faB-lV a;2+4a2 ^ ^ 

2 4 \ 7712 / 

W 9,9 ^ 1 (CT^ - 1)^ 

F2 = 27 (a/3-1) (43) 
If we remember the form of the Schrodinger-like equation which is 

+Veff^^E^, E^E^-Vo, (44) 
thus, we would write the ground-state wave-function in terms of super-potential W{x) as 



ipo{x) ^ exp{^ J W{y)dy). (45) 

We shall put the super-potential in the form of 

VF(a;) = Ci +C2tanhx (46) 

where Ci, C2 are constants, using this relation we obtain the ground-state wave-function ipQ{x) as 

Lpo{x) = e-'^^'''{coshxy^\ (47) 

There are boundary conditions as C2 > and \Ci\ < C2 such that 930(2;) — > when x — > ±00. The 
partner potentials can be given in the following manner: 

V+f{x) = W^{x) + W' (x) = Cl + Cl - [Cl - C2) sec h^x^V-2. tan hx (48) 



and 

V-j:j: (x) = (x) - W' (x) = Ci + C2 - (C2 + C2) sec h^x + V2 tan hx. 
If we show the ground state energy with Eq, we may give the expression as below 

W^{x) - W'{x) = -Visech^x + V2t&nhx - Eq. 
Now, we can match with (^0]) . one gets 

Cl + C2 = —Eo 
C2 + CI = Vi 
2C1C2 = ^2. 

Solving these equations, we obtain Ci, C2, Eq as follows 
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C2 = -(-l± 



(49) 
(50) 



(51) 
(52) 
(53) 



(54) 



and we must choose the positive sign in (IM)) because of the boundary conditions, this also leads to Vi > 0. 
The other constant Ci is given by 

C, ^ ^ (55) 



and 



En 



-1 + VI + 41/1' 

7(-l + Vl + 4Fi)' + 



(56) 



4' ' ' ' ' ' ' (-1 + VTT4VT)2' 

It is seen that two partner potentials satisfy the well-known shape invariant relationship 

V+j{x-aa)^V-jj{x-ai)+R{ai) (57) 

where oq — C2 and ai = C2 — 1. The reminder R{ai) is not depend on x and it contributes to the energy 
spectrum as 



Eq 



(58) 



fc=i 



4C| 



4(C2 - nf 



+ (C2-n)2, n = 0,l,2,... 



Eventually, using (1561) we obtain the relativistic energy spectrum for (|37p as 



^ 4(-i + iVl + 4Vi-n^' V2^ 



For real energies 1 -I- 4Vi must be positive, i.e. 



2a; + 47"^ (cr - TO2) + 



2^ , 4(a/3-l)2 



> 



(59) 



(60) 



(61) 
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and 



Vo > — ^ -2 - (l (-1 + yTTIW) - n) 



(62) 



Hereafter we shall find the wave-function (p{x). In that case, using ([59|) in (|44|) we obtain 
- If" + (V2 tanhx - (C2 + C|) sech^x) ip = ( {C2 - - — ip 

and if we use a new variable z — — tan hx in above equation and writing the function as 

.= f^V7^VVw (64) 



(63) 



^ 2 J V 2 ^ 
then we get 

(1 - z'^)P"{z) + (-2s + 2r- (2 - 2r - 2s)z)P' (z) + n{n - 2r - 2s + l)P(z) = (65) 

where 



r = 7: \ n + -{I - Vl + Wi) - ^- , (66) 

2 \^ 2^ ^ " 2 n+i(l- Vl + 4Vi) J ^ ' 

s = + ^/l + Wi) + ^ ^- \ \. (67) 

2 \^ 2^ ^ ^' 2 n+i(l- Vl + 4Vi) j ^ ' 

Thus the unnormalised wave function and upper spinor component 4>n{x) are given by 

fi±^V7i^^V7l— )(-tanM (68) 



:p 1^ tan hx [^_±}^^j P(— ) (- tan /.x) 



(69) 

where "^^^ (— tan /ix) are the Jacobi polynomials. In addition to the results here, in [2H1 H] the 

authors obtained the spectrum of the Dirac equation with scalar, vector and pseudoscalar potentials. 
Our results are consistent with |28l[T9] in case of V2 — ^ iV2. 

4.2 Example 2: non-VT symmetric vector potential 

The choice of A{x) = 6 cothcx gives a non-VT symmetric potential which is given by 

, , E ^ / c , , 2c^(toi + TO2/3) cot /i2cx \ , , 

V[x) = E tanh cx + i — - esc hex sec hex H — — — — — (70) 

V 2 — 2c(mi + m2p) + m27smh2cx/ 



where E was given in (j36p . And the mass expression reads 

M{x) — m2'y — 2m2^ csch2x. (71) 
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(72) 



Thus, A{x), V{x) and M{x) yields the effective potential given below 

Veff{x) = E'^ + sec h^cx i 1 + + (3/4+ (mi + m2/3)2) esc /i^cx 
—2cjm2(mi + m2/3) esc hex sec hex pr- + "f^m\. 

Let us take (3 as 

= (73) 
m2 

to terminate the term esc hex sec hex in (|72p . then ()72p turns into 

Veff{x) = E"^ + sech^ex ( 1 + j + sec h^cx esc h^ex — + 72m2. (74) 
To obtain a solvable effective potential, we shall add and subtract ^e^ sech^ex to (fTl)) . we obtain 

Veff{x) = E'^ (^1 - + e^ (^^ + -^-^^ sec h"^ ex +^e^ CSC h'^ ex + ^'^ ml, < a; < oo. (75) 
It is reminded that V{x) turns into 

E e 

V{x) = E tanh ex — i- esc hex sec hex. (76) 

2 

Next, we shall give the super-potential in this form 

W{x) ^ A tan hex ~ B cot hex (77) 
then we obtain the partner potentials and ground state wave-function as 

W'^{x) - W'{x) = V-j^{x) = {A-Bf + B{B- e) CSC ex- A{A + e)sech'^ ex (78) 

W'^{x) + W'{x) = V^jf{x) = {A - Bf + B{B + e) CSC ex - A{A- e) sec ex (79) 

and 

— A S. 

ipQ{x) — {cos hex) = [sin hex) <= (80) 
here ^ > and > is taken owing to the boundary conditions. Now, let us compare (1751) and (I75|) . 

iA~B)' - £;2(l-l)+7^m^ (81) 

B{B-e) = \e^ (82) 

A{A + e) = + (83) 

hence we obtain B — ^ , A = ^ — — — ^ —■ Shape invariance relation is written as 

V+j{x,ao) = V-j^{x,ai) + R{ai) (84) 
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where uq and ai are given as uq — {A, B} and ai = {A — c, B + c}. If we use the expressions E ^ E'^ — Vo, 
we find 

Err = ^("fc) ^{A-Bf -{A-B- 2cnf. (85) 
fc=i 

Finally the following relativistic energy spectrum of (j70p equals 

En = ±(5v/(7TO2)2 + (A - B)2 - (A - B - 2cn)2 (86) 

where the term inside of the square root must be positive owing to obtaining real energies. Substituting 
([55|1 in dM]) we obtain 

- ip" {x) + {{A- Bf + B{B ~ c) CSC h^cx - A{A + c) sec h^cx) ip{x) ^ {{A- Bf - {A- B - 2cnf)Lp{x) 

(87) 

and we use a new variable y = cosh2ca; and we express the function 1^9(2;) = (1 — y)^^'^{l + y) ^/'^P{y), 
then the above equation becomes 

(1 - y2)p" (y) + (-A-B-{B-A + l)y) P' (y) + n{n + B - A)P{y) = 0, (88) 

thus, wave-function is given by in terms of Jacobi Polynomials P,!^^*^ "^^^ ^^^\y) 

^„(.) = (1 - yf/^{l + y)-^/^p(^/-V2;-A/c-i/2)(y)^ (gg) 

Hence the upper component reads 

0„(a;) V"^27 - 2c(toi + mz/?) esc h2cx{l-cos h2cx)^^''{l+coa /i2cx)^'^/'=P^^/'=-i/2;-A/c-i/2) j-^^^g ^2cx). 

(90) 

Results are agree with those obtained earlier [42] . 

5 Conclusion 

In the present work, we have introduced a Hamiltonian model H which is in non-Hermitian form and 
mapped Ji into a physical Hamiltonian h. The time independent Dirac equation with effective mass in 
one dimension is related to h and transformed into the Schrodinger-like equation with the new complex 
vector potentials V{x) which are p7l) and (l76l) derived using the algebraic methods. In Ref . [I9] . the 
authors used real or pure imaginary vector potentials V{x). It is seen that composing V{x) into its 
real and imaginary components leads to more general effective potentials which are the elements of the 
Schrodinger-like equation. In example 1 and 2, terminating the imaginary part of the effective potential 
we have derived hyperbolic Rosen-Morse H-type solvable effective potential and hyperbolic generalized 
Poschl- Teller potential II potential. We note that the mass relations for each case are more general. We 
have obtained the solutions of these effective potential models using shape invariance method. We have 
seen that the real spectrum of the Hamiltonian given for solvable potentials cannot be obtained by using 
/3 = —a in Swanson Hamiltonian. Thus, the metric operator which is positive definite for the so called 
Hamiltonian can be searched in the next studies. 

We have introduced some graphs for the energy eigenvalues with respect to 777,2 • (|60p is used in figure 
1 and we note that different values of the parameters can lead to real or pure imaginary energy. For the 
red curve, the energy is real for the chosen parameters but it can be seen that between 771,2 — 4.2145 and 
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TO2 = 5.6142 we have imaginary energy values as «0. 0565786 and zO. 0310165 for the blue curve. If we 
compare these results, we see that when n takes the larger values, energy may take imaginary values for 
some specific values of TO2- When it comes to the figure 2, we have real energies for the chosen parameters 
but when n becomes larger again, the energy is imaginary for some values of TO2 which is < m2 < 1.404. 



Energy 




4 6 8 10 12 14 



Figure 1: Graph of f)6Q|l with respect to m2, for the red curve: n = 0, a = 2, a; = 3, 7 = 0.1, /5 = 6; for 
the blue curve: n = 3, a = 2, cj = 3, 7 = 0.1, ^ = 6 



Energy 




2 4 6 8 10 12 14 



Figure 2: Graph of ([86l) with respect to m2, for the red curve: n — 0,uj — 5, a — l,j — 10, 5 — 0.5; c — 3; 
for the blue curve: n = 3, = 5, a = 1, 7 = 10, 5 — 0.5, c — 3 
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